Demographic Stochasticity versus Spatial Variation in 
the Competition between Fast and Slow Dispersers 

Jack N. WaddelP, Leonard M. Sander^-^ Charles R. Doering'^''^''^''^ 

'^Department of Mathematics, University of Michigan, Ann Arbor, Michigan 48109-1043 

^Department of Physics, University of Michigan, Ann Arbor, Michigan 48109-1040 
'^Michigan Center for Theoretical Physics, University of Michigan, Ann Arbor, Michigan 

48109-1040 

'^Center for the Study of Complex Systems, University of Michigan, Ann Arbor, MI 

48109-1107 



Abstract 

Dispersal is an important strategy that allows organisms to locate and exploit 
favorable habitats. The question arises: given competition in a spatially hetero- 
geneous landscape, what is the optimal rate of dispersal? Continuous population 
models predict that a species with a lower dispersal rate always drives a com- 
peting species to extinction in the presence of spatial variation of resources. 
However, the introduction of intrinsic demographic stochasticity can reverse 
this conclusion. We present a simple model in which competition between the 
exploitation of resources and stochastic fluctuations leads to victory by either 
the faster or slower of two species depending on the environmental parameters. 
A simplified limiting case of the model, analyzed by closing the moment and 
correlation hierarchy, quantitatively predicts which species will win in the com- 
plete model under given parameters of spatial variation and average carrying 
capacity. 



Dispersal plays an important role in population dynamics models. Dispersal 
may stabilize the species' populations [TJ, allow exploration of new habitats, and 
influence community structure [2] . Despite the potential advantages of dispersal 
to a species, it may or may not be favored in a competitive environment. That 
is, environmental factors influence whether faster dispersers or slower dispersers 
have a competitive advantage which in turn influences whether dispersal rates 
evolve toward faster or slower values. For example, spatial variability tends to 
reduce dispersal rates [31 |4j [5] while temporal variability in the environment 
tends to increase dispersal rates [SI [3 [HI [H] ■ 

Demographic stochasticity is known to be important in population dynam- 
ics. In simple deterministic population models, the carrying capacity is often 
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a stable fixed point, implying populations will remain for infinite time. In 
a stochastic model, extinction, the only absorbing state, is certain at infinite 
time. Even at finite times, demographic stochasticity leads to increased extinc- 
tion risk |101 111] , particularly when populations are small or fragmented. The 
inclusion of stochasticity into non-linear mathematical models affects the mean 
dynamics, leading to, for example, generation [12] or exacerbation |13j of the 
Allee effect. It can also actualize sustained cycles in a simple predator-prey 
model (i.e. Volterra) that are not present in the deterministic model [14] . 

If demographic stochasticity can cause fundamental changes in the outcomes 
of simple population models, it is natural to wonder whether it can also cause 
significant changes in the competition of dispersers, and thus on the selection of 
dispersal rates. In this study we develop a mathematical model in which spatial 
variation and demographic stochasticity are sufficient to select either a faster or 
a slower disperser in a competitive environment, and analytically estimate the 
boundary between the two outcomes. 



1. Background 

The importance of demographic stochasticity is already evident in the clas- 
sic study of Hamilton and May who examined a discrete population mean-field 
model where all sites were equivalent and stationary, but where there was an 
explicit cost of migration in the form of a fixed probability for migrating indi- 
viduals to perish J^. In this scenario, Hamilton and May found a non-trivial 
evolutionarily stable rate of migration depending upon the probability of surviv- 
ing migration, and later with Comins [16] , upon the carrying capacity per site as 
well. Even in very poor environments, where the carrying capacity is unity and 
the probability of surviving migration is very small, the optimal migration rate 
goes to a non-zero constant, despite both spatial and temporal homogeneity of 
the environment. This was interpreted as being due to a form of kin selection, 
specifically by reduction of kin competition. 

Hastings [3], and later Dockery et al. |3], examined a deterministic, con- 
tinuous population, continuous space model of N competing species in a non- 
homogenous spatial environment modeled as 

^^^^^ = AVV(a;,t)-f M,(a;,i)(7(x)-^uj(a;,i)) (1) 

where Ui{x,t) is the population of the i-th species, Di is the dispersal rate 
of the i-th species, and j{x) is the heterogeneous growth rate. Under these 
conditions, they found that the species with the slower dispersal rate always 
drives a competing species to extinction. Dockery et al conjectured that this 
generalizes beyond pair-wise competitions, implying that dispersal rates will 
always evolve towards zero in environments with any spatial variation. 

In essence, the result of Dockery et al states that a slowly dispersing species 
exploits favorable environments without wastefuUy exploring the landscape. 
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The faster species does the complement, recklessly sending populations into 
poorer habitats. This can be considered an implicit cost of migration, rather 
than the explicit cost of migration in Hamilton and May. However, unlike 
Hamilton and May, the continuous population model disregards discreteness of 
and stochasticity in populations which drive migration rates upward. Costs of 
migration, whether explicit or implicit, drive migration rates downward. 

The Hamilton and May model (as analyzed by Comins et al) is consistent 
with the observation of Hastings and of Dockery et al in the continuum limit. As 
the carrying capacity on each site approaches infinity, the evolutionarily stable 
dispersal rate approaches zero. Although this was interpreted as the elimination 
of kin competition - if an infinite number of individuals can live on a site, then no 
one is competing for a position - it can be interpreted instead as an elimination 
of demographic stochasticity. Indeed, the results of Hamilton and May can be 
derived by considering stochasticity without reference to kin selection |17j . 

Many others have also considered the effect of stochasticity on dispersal rate 
evolution. Models studied might involve continuous [TF! or discrete (seasonal 
or generational) time dynamics UHl I^Oj , explicit migration cost [5TJ [T5J [H] , ho- 
mogeneous [50] or heterogeneous resource distribution [52], static [2D] or time- 
varying environments, stochastic environmental dynamics [191 122j. and spatial 
competition [TB] as well as resource competition. While all of these compo- 
nents are interesting, we focus here on a simpler mechanism for the selection of 
dispersal. We claim that in the presence of spatial heterogeneity in resources, 
no explicit cost of migration or temporal variation in resources is necessary for 
dispersal to be selected. Stochasticity, in essence, provides the spatiotemporal 
heterogeneity to promote dispersal, while spatial variation in resources provides 
an implicit cost of migration to dissuade from dispersal. The combination favors 
an intermediate dispersal rate. 

Cadet et al |18| examined a continuous time stochastic model with a static, 
homogeneous environment that incorporated both resource competition, in the 
form of a density dependent birth rate, and space competition, in the form of 
maximum population per site which requires newly born individuals to emigrate 
immediately. This model incorporates an explicit cost of dispersal by inclusion 
of dispersers in a pool in which they experience a higher death rate and no 
possibility of reproduction. They found that the demographic stochasticity 
promotes dispersal (distinct from spatial saturation induced emigration) when 
the patches are, on average, not filled, which allows dispersers an opportunity 
to improve their lot by leaving dense patches for sparser ones. They also found 
a complex dependence of dispersal on the dispersal pool death rate - this is 
attributed to the fact that, while a higher rate introduces more risk, it also 
increases the potential reward by reducing the saturation of patches. 

Kessler and Sander [23] investigated this phenomenon with a stochastic real- 
ization of the Hastings model. In this model, they found that the faster species 
could drive the slower one to extinction if the stochastic fluctuations were large 
(i.e., if there was a small carrying capacity) and/or the spatial variation was 
small. They offered a semi-qualtitative estimate, based on the balance between 
stochastic fluctuations and the implicit cost of dispersal, for the boundary be- 
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tween a fast species victory and a slow species victory. 

In an effort to quantitatively understand this parameter-space boundary 
between victory of a slower species and that of a faster species, we further ex- 
amine the relationship between demographic stochasticity and spatial variation 
in a pair of competing populations. We introduce, simulate, and analyze a 
model of two species which differ only in their dispersal rates while competing 
for spatially-varying resources. Wc find that, similar to the spatial model of 
Kessler and Sander, demographic stochasticity favors the faster disperser, while 
spatial variation favors the slower. 

2. Approach 

Our intent is to use a simple model which incorporates competition between 
dispersing species in a heterogeneous environment and which includes intrinsic 
demographic stochasticity. In the model, competition exists as resource compe- 
tition, modeled in the form of density dependent mortality. Two species will be 
considered at a time, each with its own dispersal rate. Though we refer to these 
rates as "faster" or "slower" , they do not refer to a movement speed, but rather 
a propensity for an individual to migrate to a new location. Spatial variation 
of resources is modeled as spatially dependent birth rates. 

Such a model must contain two environmental parameters: one which de- 
scribes the degree of spatial variation in birth rates (or, equivalently, carrying 
capacities), and one which relates to the spatial average of the carrying capac- 
ity. This latter parameter is an implicit measure of the degree of stochasticity in 
the model; as the average carrying capacity -and thus the average population- 
increases, the relative level of stochastic fluctuations decreases. 

The model has four parameters: the two dispersal rates and the two en- 
vironmental parameters. It thus makes sense to ask two questions regarding 
the competition. The first is, for two species with given dispersal rates, how 
will the result of the competition change as the environmental parameters are 
varied? That is, will the faster species drive the slower to extinction first, or 
vice versa, or will there be coexistence, or rapid mutual extinction (extinction 
being guaranteed in the long term in a stochastic model)? The second question 
is, given fixed environmental parameters, is there a best dispersal rate? 

We use simulations to suggest answers to both questions. We find that 
either the faster or the slower species will be preferred, depending upon the 
parameters larger spatial variation and/or larger carrying capacity tends to 
favor a slower disperser, while smaller values of the parameters tend to favor 
the faster disperser. Neither coexistence nor mutual rapid extinction is found. 
We also find that, for given environmental parameters, there is an evolutionarily 
stable dispersal strategy. 

Our goal is to find analytical results which explain the findings of the sim- 
ulations. This requires several compromises in the details of the model, though 
the underlying structure is maintained. Our fundamental approach is to gen- 
erate systems of diff'erential equations for the moments and correlations of the 
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populations, and then to perform stability analysis on the fixed points of the 
system. 

However, these moment equations are not closed: lower moments depend on 
higher moments, which depend on still higher moments, and so on ad infinitum. 
In order to close some of the moments and make progress, we will consider 
some simplifying limits of the model. The most dramatic of these is to assume 
extremes in the species' dispersal rates: the fast species will disperse infinitely 
quickly and the slow species will not disperse at all. This means that the fast 
individuals will sample all sites equally and will respond to the mean field. 
Importantly, this docs not mean that the variance in the fast population goes 
to zero. The zero dispersal rate for the slow species allows us to close the 
correlations between the slow population and the spatially heterogeneous growth 
rates. 

These limits leave us with only the third central moment of the slow pop- 
ulation unclosed. We show that forcefully closing this, by neglecting the third 
central moment, leads to a quantitatively correct prediction of which species is 
the victor. It also accurately predicts the trajectories over time of the moments, 
provided that the slow population does not become too small. 



3. Model 

Wc consider two species that jointly reside and compete for resources at spa- 
tially distinct sites. One species is the relatively mobile "fast" disperser, and the 
other is the relatively sedentary "slow" disperser. The space is inhomogeneous 
in the sense that the carrying capacities of the sites vary. 

The following transitions characterize the model process on L sites i, where 
1 < i < L. Fi is the number of fast-dispersers at site i, Si is the number of 
slow dispersers, 7^ is the growth rate at site i. Individiials experience density 
dependent mortality scaled by the parameter f3 = 1/n; this sets the population 
scale (which is finite), so that the carrying capacity at a site i is 7,n: 



Fi — ^ Fi + 1 Birth for fast species at site i 

Si ^ Si + 1 Birth for slow species at site i 

Fi ^ — ^ -Fi — 1 Logistic death for fast species at site i 

Si — ^ Si — 1 Logistic death for slow species at site i 

{Fi,Fj) " ~ ') {Fi - 1, Fj + 1) Dispersal for fast species from site i to site 

{Si, Sj) ^'^^^"^^ {Si - 1, Sj + 1) Dispersal for slow species from site i to site 

For example, the site-dependent growth rate 7, can chosen to be the "all or 
nothing" scenario: 

7j = 1 with probability (j) (fertile sites) (3) 
7i = with probability 1 — (j) (sterile sites) 



Site 1: Sterile site 3: Fertile 




Site 2: Fertile 



ssi 1 

Figure 1: Illustration of the system given by Eqs |2]-|3|. The large squares represent the L 
sites, and the small (filled, empty) circles represent (fast, slow) individuals. There are three 
processes: birth, death, and dispersal. Birth occurs at rate 7; = 1 on fertile sites and 7i = 
on sterile sites (i.e., birth does not occur on sterile sites). Death occurs at rate (/i + Si)/n. 
Dispersal occurs for F at rate Df and for S at rate Ds- 



The phenomena studied herein are observed to be insensitive to the specific 
choice of the ji provided that there is a tunable spatial variance (which in this 
case is controlled by the single parameter 0) . See Figure 1 for an illustration of 
this process. 

We studied this process in two ways and compare the results. First, the 
process was simulated exactly, producing realizations of the Markov process 
according to its definition. We refer to these as the direct Monte Carlo simu- 
lations (DMCS). Secondly, we analyze the differential equations of the moment 
hierarchy generated by the process. 

In order to make analytical progress we examine the simultaneous limits of 
many sites, infinitely fast dispersal for the fast species, and no dispersal for the 
slow, i.e. L — ^ oo, Df — > oo, and Ds — >■ 0. It is important to note that we 
do not take the infinite n limit in which the process would be deterministic. 
This analysis confirms that the demographic stochasticity due to finite n is 
important, as demonstrated by DMCS. 



4. Simulations 

For each particular simulation, the parameters of number of sites L, the 
fraction of fertile sites cj), the population scale n, and the dispersal rates Df 
and Ds were fixed. Initial conditions were that every fertile site started with an 
equal number of fast and slow individuals, half of the local carrying capacity for 
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Figure 2: Example mean trajectory from simulations. (/) and (s) are the mean populations 
normalized by n, and Vf and Vb are the normalized population variances, over L = 100 sites 
and 100 simulations. For these simulations, n = 50, <j) = 0.85, Df = 10, and Dg = 0. In this 
scenario, S drives F to extinction. 



each (the carrying capacity took only even values) , and that all sterile sites were 
initially empty. The time course of the populations on each site were recorded 
for each, as was the eventual winner, i.e., the species that drove the other to 
extinction (no coexistence in any realization was observed for Df ^ Dg). By 
averaging over many simulations, trajectories were generated for the moments 
of the populations (cf. Fig [2]). The fraction of victories for each species was also 
recorded for each given set of parameters. 

While varying only the fraction of fertile sites (/) and the population scale 
n, two distinct phases emerge, separated by a boundary (Figjsj). Below this 
boundary, the faster species wins, while above the slower species wins. This is 
in contrast with the results for the deterministic n — > cxi limit, where the slower 
species always wins. While DMCS recovers the deterministic result for large n, 
the faster species wins when the population scale n or spatial variance of growth 
rates Var{'-^) = — (f>) is small enough. This phenomenon was first explicitly 
noted by Kessler and Sander ^2^. 

Kessler and Sander made a semi-quantitive argument for the form of the 
fast / slow phase boundary based on the competition between the degree to which 
each population fluctuates and the degree to which each population best follows 
the contours of the food distribution. This argument determined that the phase 
boundary should scale as the following: 

Var{-f) - 1/n (4) 

where n is the average carrying capacity. The goal of the following section is to 
develop a quantitative understanding for this phase boundary. 
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Figure 3: Phase boundary between victory by the slow species (above) and victory by the 
fast species (below). The points are the result of continous-time rejection- free Monte Carlo 
simulations, with the dispersal rates Df = 10 and = 0. 



5. Analysis 

From the Master equation corresponding to ([2| and ([3| we generate differ- 
ential equations for the evolution of the moments. Let fi be Fi/n. In order to 
make the averaging precise, we consider the mean-field limit where L goes to 
infinity. We define the averages and the variances according to: 



(/)=E( lim i-X^/O (5) 
Vf - if) {f? = E( hm \ Y.{U - im (6) 



i=l 



Similar definitions hold for S and high moments of the random variables, leading 
to the following moment equations: 



dt 
dvf 
dt 



and 



{lf)-{f^)-{ff~vf (7) 
2((7/'> - (7/>(/> + (/>' + {f)vf - if) + (/>(/.) - ifs)) 
-2Df{vf ^^-^) + ((7/> + iff + vj + {fs))/n (8) 
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d{sl 

dt 
dvs 

dt 



{js)~{fs)-{.sy-v, (9) 
-2Ds{vs - ^) + {M + + v, + {fs))/n, (10) 



and the correlation equation 
d{fs 



dt 



= 2{rfs) 2{{s.f) + ifs')) + (As + Df){{f){s) (/.)). (11) 



Now we focus on the hmits of Df — )■ oo and Dg — !■ 0. In the large Df 
limit, there is rapid relaxation to Vf = {f)/n. This frees us from the need 
of a closure for the {f^) and (/^s) terms. Similarly, (/s) rapidly relaxes to 
(/s) = (/) (s) . The simulations confirm the accuracy of these approximations; 
even for Df — 10, the error was only a few percent. 

The equations for the higher correlations are readily computed, but lengthy. 
Note in particular the following: 

&l^Df{{f){s')-{s'f))+Oil). (12) 

In the large Df limit, there is also a rapid relaxation of (s'^f) to (/)(s^). This 
closure approximation was also verified by the simulations, and is found to be 
accurate even for moderate D f . Then taking -> we eliminate the Dg term 



in (101 



Next, the correlations (7s™) for integer m are calculated using the "all or 
nothing" definition of gamma ([3|. Since Ds = 0, a site with no slow individuals 
remains devoid of them for all future times. Therefore, given initial conditions 
where populations start only on the fertile sites, 

L L=cf>L L 

i—l i—1 i—(pL+l 

This applies similarly for all moments involving 7 and s. That is, (73^) = (s^), 
and so on. However, F moves about rapidly across all the sites, sampling them 
equally in the large Df limit. Thus the correlation (7/) = which the 

simulations confirm as well at even moderately large Df. 
We can then simplify the equations to the form 

^ = {<P-{.f)-{s)-l/n){f) (14) 

^ = (1 -(/)-(.) -«,/(.))(.) (15) 

^ = 2{vg + {sf + {s)vs-{s')) + {{s) + {sr+ vg + if) {s))/n (16) 
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It is now convenient to make a change of variables. Let p(t) be the mean 
of the slow population on the fertile sites. Since the slow population is zero 
on the sterile sites, (s) = (fjp. The variance Vs is comprised of two components 
- the variance on the fertile sites Vp, and the variance between the populated 
fertile sites and the abandoned sterile sites. Thus, = (j)Vp +p'^(j){l — <j)). For 
convenience, we apply the same scaling with the fast population so that (/) = cfyq 
defines the new variable q{t). 

Let ^ be the 3rd central moment of the slow population on the fertile sites. 
That is, ^ = {{pi — p)^), where pi is the random value of the slow species 
population on each fertile site, so that p, defined above, can be written as 

P = (Pi)- 

The equations then become 



dq 

dt 
dp 

di 
dvp 
~dt 



<t>{l-q-p-\l{n<t>))q (17) 
{l-(f)q-p-Vp/p)p (18) 
2{l-(f)q-2p+ l/{2n))vp + {1 + p + (l)q)p/n - 2^ (19) 



The remaining challenge is the third moment ^, which is not closed. In 
circumstances when p is large, such as when the slow population is dominating, 
we expect that ^ is small compared to the other terms and can be safely dropped. 



The simulations bear out this expectation. Dropping the last term in ( 19 ) leads 
to an accurate prediction when the slow population wins (see Figure |4|). When 
the slow population is near extinction, ^ is not small compared to the other 
terms in the Vp equation and it cannot be safely neglected. Neglecting ^ in 
these situations leads to instability of the ODE system: the evolution drives Vp 
to unphysical negative values and p diverges (see Figure [5]). 

Thus, while neglecting ^ does not always reproduce the results of the stochas- 
tic system, it still leads to an accurate prediction of the final outcome. When the 
ODEs predict that the slow population will win, then the ODEs are well behaved 
and produce good quantitative approximations. When the ODEs demonstrate 
instability, then ^ is not small, indicating that the slow population is going ex- 
tinct, i.e. that the fast population is winning. That is to say, the model always 
offers a quantitative prediction of the binary variable describing whether F or 
S is the ultimate victor, even if it does not always quantitatively predict the 
moments as a function of time. 



When is neglected to close the system, Eqs. (17 - 19) have four pertinent 
fixed points. A finite stochastic system like this will eventually end in total 
extinction, the only absorbing state, at a time that is exponentially large in 
the population. Rate equations are not sensitive to these large deviations, and 
stochastic differential equations get quantitatively wrong results, so more so- 
phisticated techniques would have to be used to determine the extinction time 
|24j. This ultimate extinction occurs in the simulations, which of course require 
finite L, but because L is large the extinction time is very long compared to 
the time scale of a single-species victory. The ODE analysis is for the case 
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Figure 4: Example of the mean populations from the ODEs and the simulations when S drives 
F to extinction. {/} and (s) are the mean populations, and Vf and Vs are the population 
variances. The solid lines are given by solving the ODEs (or, in the case of Vf, by the closure 
used in the ODEs), and the points (circles for F and crosses for S) represent the moments over 
100 simulations. For these simulations, n = 40, L = 100, (f) = 0.5, Df = 10, and Ds = 0.001. 
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Figure 5: Example of the mean populations from the ODEs and the simulations when F drives 
S to extinction Note that the ODE trajectory initially matches the simulations until suddenly 
it diverges and S becomes negative, at which times the ODE integration is ceased. Although 
the trajectories diverge at this point, the ODEs have made a prediction for the outcome: S 
will go extinct. For these simulations, n = 40, L = 100, (j) = 0.99, Df = 10, and Ds = 0.001. 
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Figure 6: Fixed points of the rate equations versus the quasi-steady state moments from the 
simulations. In each case, the competing species population was set to zero, so the equation 
follow the competition-free trajectory to the quasi-steady state. The mean and variance 
of both the fast and slow species' populations were well matched by the rate equations, 
particularly as the population scale n increased. These were computed with Df = 10, Ds = 0, 
L = 100 and </< = 0.95, for n ranging from 10 to 160. 



that L — > cx), in which case the extinction time diverges, so this aspect is not 
contained in the ODE modeL 
The four fixed points are: 

(.<l,P,Vp)= (0,0,0) Extinction (20) 

(1-^,0,0) Victory by fast (21) 

= (0,i(3-f y^),|(l-y^+^)) Victory by slow (22) 
= (0^(3 1(1 + V^+t)) Saddle (23) 

These expressions are derived for large n, and only make sense for n > 8. Neither 
the extinction nor fast-victory fixed points are ever stable since both have two 
strictly positive eigenvalues. 

The fast species achieves victory with some parameters in the stochastic 
system, but because this requires a small slow population and ^ can no longer 
be safely neglected, the rate equations do not "find" this result. Despite the 
inaccuracy in the prediction of the stability, the fixed points do give the quasi- 
steady state value achieved in the simulations (to within an error of 1/L. See 
Figure [6]). 
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Figure 7: Phase boundary between victory by the slow species (above) and victory by the 
fast species (below). The points are the result of continous-time rejection- free Monte Carlo 
simulations with dispersal rates Df = 10 and Ds = 0. For each set of parameters, 100 
simulations were conducted. The points indicate the parameters which generate a 50/50 split 
for victories between S and F, and the error bars represent a 95% confidence interval. The 
solid line is the result of integrating the differential Eq. \19\ , and then decrementing n for 
each value of . The dotted line marks the stability threshold of the s-victory fixed point l |24[ | 
of the differential equations, which is a lower bound for the phase boundary. The differential 
equations match the simulations well, particularly in the high carrying capacity limit. 



The slow-victory fixed point has two eigenvalues that are always negative. 
The remaining (largest) eigenvalue is also negative when 

2 2 

" ^ 0(1 - 0) " V^)' ^^^^ 

equivalent to the expression determined in the Kessler and Sander model. Again, 
the value of the stochastic quasi-steady state is achieved by the equations. 

The final fixed point, with at least one positive and one negative eigenvalue, 
serves to split the phase space into the domain of attraction for the slow- victory 
fixed point. It is important to note that when the slow- victory fixed point 
is stable, it is the initial conditions combined with the extent of the basin of 
attraction that determine whether a specific trajectory will result in a slow or 
fast victory. That is, n > is a necessary, but not sufficient, condition for 

the slow species victory. Thus the stability threshold n = gives a lower 

bound to the phase boundary (see Figure [7| . 



It is noteworthy that the stability condition ( 24 1 yields a relationship similar 
to that estimated by Kessler and Sander for the boundary. The estimate goes 
as follows: by taking the difference in the rate equations for the means, we find: 
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= (/ - s> [(7> - (/ + .'^>] + [Cov(7, /) - Cov(7, s)] - [vf - v^] (25) 

where Cov(7, /) (7/) — (7)(/). If the two populations are the same, the sign 
of the change depends on the competition between two differences. The fast 
population smoothes its profile, reducing both its variance and its covariance 
with the 7i. Approximating these as zero leaves the difference between Vs and 
covariance of s and 7 to determine the advantage. The population s is of order 
7i, so Cov(7, s) Var(7). When is near unity, the variance Vs ~ 4>/n due 
to demographic stochasticity. Equating Cov(7, s) and Vs implies Var(7) = 
which in turn implies; 

n = (/)/Var(7) (26) 



which is similar to the stability condition given by ( 24 1 when is near unity. 



6. Summary and Discussion 

We have developed a model to investigate the effect of intrinsic demographic 
stochasticity on competition between dispersers in a heterogeneous environ- 
ment. In order to get analytic results, we examined extreme limits in dispersal 
rates, which allowed the closure of most of moment equations, leaving only the 
third central moment of S unclosed. Simply neglecting this moment led to ac- 
curate quantitative predictions of the population trajectories over time when 
S remained large, and regardless of the size of S, to an accurate quantitative 
prediction of the outcome, i.e. which species drives the other to extinction. 

Via simulation and analysis we have confirmed the observation that intrin- 
sic demographic stochasticity can favor the fast dispersal in head-to-head com- 
petition with a slower competitor. This can be intuitively interpreted in the 
following way: demographic stochasticity leads to sites which are, by chance, 
inhabited by fewer individuals than the carrying capacity allows. A fast dis- 
perser is more likely than a slow disperser to happen upon these sites, after 
which it can take advantage of the remaining capacity |21j . 

This is to be contrasted with the deterministic results of Hastings and of 
Dockery et al. that spatial variations in fecundity drive dispersal rates lower 
[21 H] due to an implicit cost of dispersing for the faster disperser when it leaves 
favorable habitats. In these deterministic works, dispersing individuals have 
no opportunity happen upon sites that are by chance favorable. The cost of 
dispersal without benefit drives dispersal rates to zero. 

Explicit dispersal costs also drive down dispersal rates [HI [16] . As dispersal 
benefits are removed, for example by reducing kin competition or demographic 
stochasticity by increasing the carrying capacity, slower dispersal is again fa- 
vored. If dispersal benefits are held constant as dispersal costs are reduced, 
either by lowering the explicit dispersal cost or by reducing the spatial variation 
in the environment, higher dispersal rates are favored. 
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Other factors are known to influence dispersal rates. Kin competition [T5] . 
habitat extirpation risk [TCI and temporal habitat variabihty O [3 [H E] favor 
faster dispersal. 

In the absence of factors detrimental to dispersal, dispersal rates in a stochas- 
tic model increase. In our model, if the environment is spatially homogeneous, 
dispersal is not a neutral factor, as it was in Hastings model. There is no cost 
to dispersal, implicit or explicit, but there is a benefit to dispersal. This leads 
to a faster disperser always defeating a slower disperser (data not shown; for 
finite L, the two species reach parity asymptotically if both dispersal rates are 
much faster than the birth rate, though cohabitation is not observed). 

When dispersal-promoting and dispersal-suppressing factors compete it is 
expected that an optimal dispersal rate might exist, the rate that best balances 
the factors and thus defeats any competitor with either a higher or a lower 
dispersal rate. There is numerical evidence in these simulations (away from 
the extreme limits in the analysis section) that this is the case in the model 
introduced here as well. The boundary, which has the same qualitative form for 
dispersal rates less extreme than the limits taken in the analysis section (Figure 
[8| , implies that the optimum dispersal rate decreases when the population scale 
n or the variance of the site quality Var(7) increases, since the slower species 
wins in these limits. We see this explicitly by varying both dispersal rates, and 
discovering that, for the parameters used, the species with D « 0.35 can both 
invade a resident population and resist invasion by a mutant (Figure [9|. 

While it is challenging to make quantitative comparisons between the pa- 
rameters in our model and empirical observations, in some cases the predictions 
seem to agree with real data. For example, Lowe |25j compared the long distance 
dispersal rates of salamanders measured by kurtosis of movement distribution 
with the goodness of habitat measured by body condition of the salamanders. 
It was found that dispersal decreased when the average body condition was high 
(analogous to large n) or when body condition variance was high (analogous to 
large Var(7)), precisely as expected. 

It is also of interest to note that the system studied in this paper is one 
for which the effects of stochastic fiuctuations have a qualitative impact on the 
result for arbitrarily large population scale n. While it is true for fixed spatial 
variation Var(7), n may be chosen to be large enough to recover the deterministic 
result (viz., that the slow species wins), there is no value of n which is large 
enough to always ensure a deterministic result for arbitrary Var(7). That is, 
for any n, there exists a Var(7) < 2/n small enough such that the effects of 
demographic stochasticity leads to a fast species victory. 
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Figure 8: Phase boundary between victory by the slow species (above) and victory by the 
fast species (below) for three pairs of values for dispersal rates Dr and Ds ■ The qualitative 
forms of the boundaries are similar, and also similar to figure Isl However, the boundary 
shifts as the dispersal rates Dj and Ds are varied. This is to be expected if there exists an 
optimal dispersal rate which depends on the population scale and spatial variation in habitat 
quality. For example, if one species has a dispersal rate that is fixed at the optimum for a 
given population scale and spatial variation, then that point on the phase diagram should 
be in the slow-victory phase if the competing species has a faster dispersal rate, but in the 
fast- victory phase if the competing species has a slower dispersal rate. 



16 




ss 



ai 0^ 0.3 04 OA ae o.t o.s aa i 

Pmldenl D 



Figure 9: Pairwise invasibility plot of dispersal. The resident species begins at carrying 
capacity on each fertile site, and the invading species begins with one individual on each 
fertile site. Color represents the fraction of trials in which the invading species goes extinct. 
An optimum dispersal of D* « 0.4 is evident. In this example, 100 simulations were run for 
each, point, and n = W, L = 500, and <j> = 0.5. 
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